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We propose to use the spin-blockade regime in double quantum dots to reduce nuclear spin polar- 
ization fluctuations in analogy with optical Doppler cooling. The Overhauser shift brings electron 
levels in and out of resonance, creating feedback to suppress fluctuations. Coupling to the disordered 
nuclear spin background is a major source of noise and dephasing in electron spin measurements in 
such systems. Estimates indicate that a better than 10-fold reduction of fluctuations is possible. 



A common goal across all areas of experimental physics 
is to reduce noise and fluctuations. While this can often 
be achieved by cooling the entire system, in many situ- 
ations only a few relevant degrees of freedom need to be 
"cold." Directed cooling methods have been developed 
to specifically target these experimentally important de- 
grees of freedom, allowing the relevant subsystem's tem- 
perature to be reduced by many orders of magnitude with 
little effect on the rest of the system. One very impor- 
tant example of directed cooling that has led to great 
successes in the field of atomic physics [1] is the method 
of Doppler cooling [2] , used to cool the velocity distribu- 
tion of atomic and molecular gases. Cooling is achieved 
by absorption of near-resonant monochromatic radiation, 
with the velocity-dependent Doppler shift enhancing ab- 
sorption for faster particles and suppressing it for slower 
particles. In solid state systems, laser-like cooling was 
proposed for nanomechanical resonators [3^, 0, [H , and re- 
cently demonstrated for a superconducting qubit . 

Advances in semiconductor technology have opened up 
many new possibilities for controlling microscopic degrees 
of freedom in solid state systems. Experiments in the so- 
called spin-blockade regime of double quantum dots 
have demonstrated control over the charge and spin de- 
grees of freedom of individual electrons, as well as con- 
trol over the spin degrees of freedom of lattice nuclei due 
to their hyperfine coupling to electron spins This 
same coupling of electron and nuclear spins also gives rise 
to the primary source of noise and dephasing in electron 
spin measurements due to the disordered nuclear spin 
background 13, As a result, there is now a great im- 
petus to better understand how to control and to cool the 
nuclear spin degrees of freedom in such systems H, [l^ . 

The key idea behind Doppler cooling is that the pro- 
cess that leads to cooling (absorption) receives feedback 
from the degree of freedom that is being cooled (veloc- 
ity). Here we propose an analogous scheme for cooling 
the nuclear spin polarization degree of freedom in a dou- 
ble quantum dot. Using a simple model of spin-blockaded 
electron transport, we identify a window of detuning 
where nuclear spin evolution is controlled by the balance 
between competing rates of nearly resonant spin flip-flop 
processes in opposite directions. These rates are sensitive 
to polarization via the Overhauser shift, which brings 
electron singlet and triplet states in-to and out-of reso- 




FIG. 1: Resonant suppression of nuclear spin fluctuations. 
Net spin flip rates are peaked near the resonance of the elec- 
tron triplet states with the broadened singlets | S ) and \S' ). 
Dashed lines indicate the energies of the electron triplet levels 
\T±) in the steady state. A fluctuation of polarization that 
changes the triplet level splitting leads to an imbalance of 
rates that drives polarization back to the steady state value. 
Fluctuations of either sign are suppressed (see times ti,2)- 

nance. The interplay between opposing processes gives 
rise to a rich phase diagram that exhibits hysteretic be- 
havior. In direct analogy with Doppler cooling, we show 
that nuclear spin fluctuations are strongly suppressed in 
phases dominated by nearly resonant spin flips. We quan- 
titatively analyze this cooling effect and discuss possible 
ways to optimize and employ it in experiments. 

The key physical principle responsible for cooling in 
the spin-blockaded double dot is illustrated in FigiH The 
two-electron triplet states | T± ) are split by the Zeeman 
energy Ez = :^giJiB[B + as{t)]^ where B is the applied 
magnetic field, and as{t) is the Overhauser field felt by 
electrons due to hyperfine coupling to the nuclear spin 
system with instantaneous mean polarization s(t). The 
electron g-factor and hyperfine coupling in GaAs quan- 
tum dots are ^caAs ^ —0.4 and aN — 5 T where N is 
the number of nuclear spins in the dot. Because a < 0, 
nuclear spins polarized against the external field increase 
the total field seen by the electrons. 

Transitions between electron triplet and singlet states 
occur due to the hyperfine interaction. Each transition 
from I T_ ) to | S*' ) is assisted by the flip of one nuclear 
spin from up to down. Similarly, transitions from | ) to 
I S ) are assisted by spin flips from down to up. The rates 
of these transitions are greatly enhanced when electron 
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singlet and triplet levels approach resonance. 

In the steady state, the total rates of nuclear spin flips 
from down to up and from up to down are equal. When 
this steady state is perturbed by a fluctuation of nuclear 
spin polarization in the direction opposite to the exter- 
nal field, electron Zeeman energy increases. For electron 
states detuned as in FiglU this brings the levels | T+ ) 
and I S ) closer to resonance and pushes | T_ ) and | S' ) 
further out of resonance; the rate of nuclear spin flips 
into the direction of the external field becomes greater 
than the rate of opposing spin flips and polarization is 
driven back toward the steady state value. Fluctuations 
in the opposite direction are suppressed analogously. For 
materials where a is positive, cooling will be achieved for 
energy detuning on the opposite side of the resonances 
where the sign of feedback is reversed. 

To understand this cooling effect quantitatively and to 
determine its experimental signatures, we now present 
a detailed analysis. We start by generalizing the model 
put forth in Refs.0, [1,1,^0] to describe spin blockade of 
their double quantum dot devices in terms of two-electron 
states (0,2)s, (1, l)s and (1, The states {nL^nR)s/t 
are labeled by electron occupation numbers ul and ur 
of the left and right dots. Subscripts denote the two- 
electron spin state, which can be singlet or triplet. 

Recent experiments in lateral double dots [§, [l^ were 
made in the regime where the singlet states (0, 2)^ and 
(1, l)s can be tuned through an avoided crossing. Orbital 
hybridization makes the labeling scheme based on dot 
occupancy not suitable for the singlets. Instead we will 
use the labels S and S' for singlets, and To,± for triplets 
as shown in FiglU Throughout this work we will take[16i] 



Es'^s = U^V^T^ - A), Et, = 0, Et^ = ±Ez^ (1) 

where t is the tunnel splitting, and A is the detuning 
from the singlet energy crossing. 

Transport occurs as a sequence of tunneling events. 
Charge moves from the source to the left dot, through the 
right dot, and finally to the drain. Starting with one elec- 
tron localized on the right dot, the second electron tun- 
nels into the left dot to form one of the five two-electron 
states {S^ S' ^To,T±}. The probabilities for filling these 
states are determined by the relative rates of tunneling in 
from the source lead. Because the three triplets have the 
same (1,1) orbital wavefunctions, the rates for tunneling 
into any of these states are the same. 

The mechanism of charge transfer to the drain lead is 
very different for the singlet and triplet states [7] , as the 
coupling of the (0, 2) state to the drain is much stronger 
than that of the (1,1) states. Charge moves easily to 
the drain from the states S or S' ^ since they both have 
a finite component in the (0, 2)^ state. When one of the 
triplet states is filled, however, the out-flow of charge is 
blocked; Pauli's Exclusion Principle and the high orbital 
excitation energy prevent the system from forming the 



(0, 2)t state. Residual current is due to hyperfine fiip-fiop 
scattering into the singlet states S / S' , exchange with the 
leads, or tunneling through virtual excited states. 

The rate of spin-flips due to hyperfine coupling is gov- 
erned by the near resonance of the triplet states T± 
with the singlet states S and S' . We describe the nu- 
clear subsystem by the populations N± of the up and 
down nuclear spin states, neglecting their spatial varia- 
tion. The energy-dependent hyperfine spin flip transition 
rates W^^ are the sum of rates for transitions to S and 
S' ^ calculated using Fermi's Golden Rule: 
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where A is the orbital overlap factor between the triplet 
state T± and the singlet £± is the energy difference 
between T± and 5, and 7 is the width of S. We as- 
sume a Lorentzian lineshape to allow explicit calculation. 
Primed terms represent the same quantities, calculated 
with respect to singlet state S' . Throughout this paper 
we take A = A' = 1 without loss of generality. 

The net spin flip rates are controlled by the compe- 
tition between the rates of hyperfine spin flips and non- 
spin- flip tunneling processes in each channel [14]. Assum- 
ing a single energy independent non-spin-flip escape rate 
W^"^ and equal tunneling-in probabilities of 1/4 for each 
of the triplets, the net nuclear spin flip rates are given by 



r± = 
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where Iq is the inverse of the average dwell time in To,± 
with To"^ = W'"" and = {W''' + W^^). We assume 
that the time scales for loading and for decay of the 
singlets are much shorter than the average triplet life- 
time, and that the To state does not decay by hyperfine 
transitions due to lack of energy conservation. The spin 
flip rates r± receive feedback from nuclear polarization 
through the combination of the resonant energy depen- 
dence (|2]) and the Overhauser shift of electron levels [jj . 

We now study the flow of the average nuclear polar- 
ization s = N-^ — N- . Because every spin flip changes s 
by two, polarization evolves according to 



F(s) = 2(r+-r_)-rreis, 



(4) 



where Frei is a phenomenological relaxation rate due to 
nuclear spin diffusion and thermalization. In the absence 
of a non-spin-flip escape path, VK^" = and every spin 
that enters the dot must flip its spin; r+ = r_ and no 
spin can be pumped. When W^^ is finite, a net polar- 
ization rate V{s) 7^ is achieved by tuning near one 
of the singlet-triplet resonances where spin flip processes 
are favored in one channel while non-spin-flip tunneling 
is favored in the other. Polarization flows toward fixed 
points V{s^) = which give the steady state values of 
nuclear polarization. 
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FIG. 2: Steady state polarization = versus magnetic 

field found numerically using dH), (|3]), and {|2]) with N ^ 10^, 
t = 35/ieV, A = 3/ieV, W'^ = 10S"\ and Trei = 0.5s"\ 
Blue circles denote stable fixed points V\s*\ B) < 0, while red 
crosses denote unstable fixed points V\s*; B) > 0. Phases of 
type I are relaxation dominated, while phases of type II are 
resonance dominated. Left Inset: Width of the nuclear spin 
distribution relative to the thermal width aTh = a/^V. Right 
Inset: Hysteretic current with magnetic field sweep, calcu- 
lated using lo defined in (|3| (see discussion before Eq.JS))). 
Arrows indicate sweep direction. 



With the help of the stability condition, V'(s*) < 0, 
we find that V{s) contains 1, 2, or 3 stable fixed points 
depending on the values of parameters such as detuning, 
magnetic field, and relaxation rate. In Fig|2]we plot 5* 
as a function of magnetic field at constant detuning with 
realistic parameter values [13] • 

Examining FigO we identify two distinct types of 
phases, unpolarized (labeled I) and polarized (labeled 
II). Phases of type I are relaxation dominated, and are 
characterized by nearly zero average polarization. Spin 
flip transitions are off-resonant in these phases, and nu- 
clear spin dynamics are dominated by noise and relax- 
ation. Phases of type II are resonance or polarization 
dominated. In these phases, spin flip transitions are res- 
onant in both directions, allowing for the build up of nu- 
clear spin polarization and resonant suppression of fluc- 
tuations. The existence of such phases is made possible 
by the close proximity of both resonances S and S' . 

With moderate relaxation Frei, we find that the un- 
polarized state is stable at zero applied magnetic field 
(phase I). This phase remains stable as magnetic field is 
increased, until a bifurcation point is reached (approxi- 
mately B = 0.5 T for the case shown). At this point, 
the Zeeman-split triplet levels reach far enough into the 
tails of the singlet resonances S and for the hyperfine 
spin-flips to overcome relaxation and trigger an instabil- 
ity towards a spin-polarized state (phase II). 

After entering phase II, polarization changes smoothly 



from negative to positive values as B is increased. At even 
higher fields, a second bifurcation point is reached where 
the polarized phase terminates (5 ^ 1.2 T in the exam- 
ple shown). Polarization flows back to and the system 
returns to the relaxation-dominated phase V . On sweep- 
ing back towards lower fields, the system follows a differ- 
ent path. This hysteresis is manifested in a double step 
of current for the magnetic field sweep described above 
(see Fig 12] inset), which is in striking resemblance to the 
current instabilities observed in Ref.[7]. 

To estimate nuclear spin fluctuations in the phases 
I and II, we shall analyze the behavior of the proba- 
bility distribution p{s^t). Sequential electron transport 
through the spin-blockaded double dot is accompanied 
by a stochastic series of nuclear spins flips that change 
the net polarization in steps of ±2. The polarization 
variable s executes a directed random walk on the in- 
terval [-N,N], where N ^ 10^ - 10^. We describe this 
stochastic dynamics by the Fokker-Planck equation : 

where V{s)At = (As) and D{s)At = l/2{As'^), with 
V{s) = 2 (r+ - r_), and D{s) = 2 (r+ + r_). Here At 
is a short interval on the order of one electron transit 
time. 

In the simplest model of relaxation, polarization decays 
as the result of random flipping noise: with microscopic 
rate Frei, a random spin is selected and its z-component 
flipped. The spin-flip rates ([3]) are then changed to 



f ± = r± + TreiiV^/iV, N± = ^{N±s). 



(6) 



This model reproduces the phenomenological relaxation 
term — Freis used in Eq.(|4j) and accounts for broadening 
of the distribution due to coupling to the thermal bath. 

The steady state of ([5]) is of the form 
exjp J^^V{s')/ D{s')ds' . One can easily check that 
this distribution is peaked around the stable fixed points 
V{s*) = 0, V^{s*) < found earlier. When fluctuations 
of s about s* are small (i.e. away from bifurcation 
points), we linearize V/D and obtain a Gaussian steady 
state distribution peaked at 5 = 5*, with rms width 
(TP = ^/D{s*)/V'is*). 

In the absence of external driving, the thermal noise 
and nuclear spin diffusion responsible for Frei cause the 
system to relax to the fully random high temperature 
distribution with zero mean {s) = and rms width 
(jTh = y/ (5^) = VN [13. As current is restored and 
as the spin flip processes are brought closer into reso- 
nance, the net suppression of fluctuations is determined 
by a competition between confinement due to the reso- 
nant spin flip dynamics and broadening due to Frei- 

The ratio dp /aTh serves as a quantitative measure of 
the strength of cooling. As Fig|2] illustrates, phases of 
type I are characterized by dp /aTh ^ 1 and thus do 
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FIG. 3: Polarization steady states (a) and widths (b) vs. de- 
tuning. Notation is the same as in Fig|2] Coohng is most effi- 
cient near A 7, where both resonances strongly contribute 
to the suppression of ffuctuations. Parameters: N = 10^, 
W'"" = 10^ s-\t = 50/ieV, 7 = 7' = 0.1 /ieV, and B = 0.5 T. 
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FIG. 4: Current /(s), Eq.Q, as a diagnostic of nuclear polar- 
ization. Results for P — are shown, when the non-spin-ffip 
escape rate W^^ is due to tunneling back to source lead. Note 
the non-monotonic dependence of current on polarization. 
Parameters: N 10^, t 24/ieV, A 4/ieV, B 0.2 T, 
W'"" = 10^s"\ and Frei = 0.1 s"\ 



not exhibit cooling. The polarized phase II, however, 
displays a significant suppression of fluctuations. Away 
from the edge of the instability, the strength of cooling 
is insensitive to magnetic field. 

In contrast, the efficiency of cooling depends very 
strongly on detuning. In Figl3] we plot polarization and 
cTp /(jTh as a function of detuning in each of the three 
phases I, II, and II' present for the chosen parameters. 
Maximal cooling occurs at small positive detuning 



7, 
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where the system can benefit from the resonant character 
of spin flip rates for both the T+ and T_ states simultane- 
ously (see Figs. [T] and [8)3). With both resonances 5* and 
contributing to the Fokker-Planck effective restoring 
potential, a greater than 10-fold suppression of thermal 
fluctuations can be achieved (see FiglSb inset). 

How are these phenomena reflected in device cur- 
rent? Depending on the mechanism(s) responsible for 
the non-spin- flip escape rate VK^", the current / trans- 
mitted through the device may differ from the rate /q in 
(|3|). Only electrons that move from source to drain count 
toward the flow of current; those that jump back to the 
source do not contribute. Thus we introduce a parame- 



ter < /3 < 1 that specifies the fraction of W^^ due to 
non-spin-flip escape to the drain, yielding 
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When = 1, all electrons exit to the drain and I = Iq. 
We stress that the spin flip dynamics depend only on Iq 
and not on the mechanisms responsible for W^^ or on /3. 
Figure 2] illustrates that current ([8]) has a complicated 
non-monotonic dependence on 5, but nonetheless can be 
used as an indicator of polarization. 

We have shown that resonant hyperfine spin flips can 
be used to significantly suppress nuclear spin fluctuations 
in spin-blockaded double quantum dots by a mechanism 
analogous to Doppler cooling. Optimal cooling occurs 
at small positive detuning A* ^ 7 where both triplet- 
singlet resonances provide strong feedback through the 
Overhauser shift. The suppression of nuclear spin fluc- 
tuations will manifest itself dramatically in electron spin 
resonance (ESR) measurements [11]. ESR experiments 
will both benefit from cooling through increased dephas- 
ing times and reduced fluctuations in Zeeman energy, and 
serve as a powerful way to detect it. 
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[16] In previous work [14] we discussed transport in the regime 
of large detuning, where the energy separation of the 
(0, 2)s and (1, l)t states is much larger than that of the 
(1, l)s and (1, l)t states. 
[17] This result is true for ID (Ising) spins. The width of the 
disordered state of 3D spins is (Jth = a/ N/3. 



